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Abstract 

For the interactions between a He(2 3 S') atom and a He(2 3 P) atom for like isotopes, we report 
perturbation theoretic calculations using accurate variational wave functions in Hylleraas coordi- 
nates of the coefficients determining the potential energies at large internuclear separations. We 
evaluate the coefficient C3 of the first order resonant dipole-dipole energy and the van der Waals 
coefficients Cq, C$, and C\$ for the second order energies arising from the mutual perturbations 
of instantaneous electric dipole, quadrupole, and octupole interactions. We also evaluate the co- 
efficient C9 of the leading contribution to the third order energy. We establish definitive values 
including treatment of the finite nuclear mass for the 3 He(2 3 5)- 3 He(2 3 P) and 4 He(2 3 5)- 4 He(2 3 P) 
interactions. 

PACS numbers: 34.20.Cf,31.15.Ar,31.15.Ct 
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Recently, there has been considerable interest in the study of helium dimers and ultracold 
helium collisions associated with metastable helium atoms 

Molecular lines of the helium dimers associated with the He(2 3 S')-He(2 3 P) asymptotes have 
been produced by photoassociation of spin-polarized metastable helium atoms He(2 3 Si). 
In the purely long-range 0+ potential well (with a minimum internuclear distance reaching 
150 bohr) associated with the He(2 3 Si)-He(2 3 P ) asymptote, bound states of the helium 
dimer can be used in measurements of the s-wave scattering length for collisions of 

two 4 He(2 3 Si) atoms 0, Q|. Molecular lines to the red of the D 2 atomic transition in helium 
dimers associated with the He(2 3 S')-He(2 3 P2) asymptote are of interest for control of the 
scattering length using an "optical Feshbach resonance" 0, E|. Published data are scarce 
on the long-range part of the He(2 3 S')-He(2 3 P) potential energies, which determine the 
energy level structures of the ultralong-range dimers formed in photoassociation of ultracold 
metastable helium atoms. 

In this paper, we present perturbation theoretic calculations of the coefficients determin- 
ing the potential energies at large internuclear separations using accurate variational wave 
functions in Hylleraas coordinates. We evaluate the coefficient C3 of the first order resonant 
dipole-dipole energy, and the van der Waals coefficients C 6 , C 8 , and C\q of the second order 
energies arising from the mutual perturbations of instantaneous electric dipole, quadrupole, 
and octupole interactions. We also evaluate for the third order energy the coefficient Cg 
of the leading term. We establish definitive values including treatment of the finite nuclear 
mass for the 3 He(2 3 S')- 3 He(2 3 P) and 4 He(2 3 S')- 4 He(2 3 P) interactions. Definitive values for 
the C*6, Cg, Cg, and C10 coefficients are established. 

In this work, atomic units are used throughout. For the He(2 3 S')-He(2 3 P) system, the 
zeroth-order wave function appropriate for the perturbation calculation of the long-range 
interaction can be written in the form [12J: 

* (0) (M, ±) = -±=[* n >)*n b (lM; p) ± MUp)*n 6 (lM; <x)] , (1) 

where ^ Ua is the metastable He(2 3 S) wave function, ty Ub is the wave function for the He(2 3 P) 
atom with magnetic quantum number M, a and p are the collection of coordinates of the 
two helium atoms in the laboratory reference frame, and the " ± " sign indicates the gerade 
and ungerade states. The corresponding zeroth-order energy is = E^> + E^J . At 

large internuclear distances, the interaction potential between the two helium atoms can be 
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expanded as an infinite series in powers of 1/R 12 



v = EE 



v, 
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where 
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Tfp(a) and T_J(p) are the atomic multipole tensor operators defined as 
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and 
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where Qi and qj are the charges on particles i and j. The coefficient i^^ L is 
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+ L\ /£ + L\l 1/2 
' + //J \L + n)_ 

and (£, Lj. . .) = (2£ + 1)(2L + 1) . . . . According to the perturbation theory, the first-order 
energy [12j can be written as 

C 3 (M,±) 



i? 3 



where 
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c s (o,±) = ±Yl(^ a (o-)iiEQ^(^)ll^(i;o-))r 
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and 



Air 

C 3 (±l,±) = Ty|(*n a (o-)||X:^n(^)||^(l;o-))r- 
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The second-order energy is 



^ (2) = -£' E E 



I (jg]V]Xn. (L S M S ; cr)uj m (L t M t ; p)) p 



(10) 



n s n.t L S M S LtMt 

where Xn s (-^s^f s ; ff)u nt (L t M t ; p) is one of the intermediate states with the energy eigenvalue 
-E„ ant = E^J + Ej^>, and the prime in the summation indicates that the terms with E naUt = 
E^J should be excluded. Substituting (J2J into (fTUJ) . we obtain 

B x ±B 2 



= -£' E E 



n s r»t L a M s L t Af t -^n s n t — ^n, 



(0) 
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with 

Bi = |(v&„ a (o-)vI/„ 6 (lM;p)|y| Xns (L s M s ;o-V nt (LiMi;p))| 2 , (12) 

and 

B 2 = {^ n S(r)^ n SlM ] p)\V\xnXLsM s - 1 a)u nt {L t M t]P )) 

x (* no (p)* nfc (lM; <r)\V\ X n s (LsM s ; cr)uj nt (L t M t ; p)) . (13) 

B 2 in Eq. (jl3|) is the exchange interaction of the two states He(2 3 S) and He(2 3 P). After 
applying the Wigner-Eckart theorem, we have 



Ei sr^ \p -Si __ \p Ci(L, L', L s , L t , M) 

2^ 2^ „ F (0) ~~ 2^ T>2L s +L+L'+2 ' v ' 

L S M S L t M t -t^ns-nt ~ -^riant LL'L s L t 

with 



C x {L,L',L.,Lt,M) = ^G , 1 (L,L' J L s ,L t7 M)F^(L,L' J L s ,L t ). (15) 

In Eq. (|15jh is the angular-momentum part and is the oscillator strength part. Their 
expressions are 

G[(L,L',L S ,L U M) = (-1) L + L '(L,L')V> £ 

M s A/t 

x( 1 L L ')( 1 L ' («» 
V -M M s MJ V — M M s M, / 



and 



with 
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(vl> nQ (Li; <r) 1 1 £ llXn. (£.; »)> 

(*^(L 2 ;ff)||$3Q < (7f^(ff i )||xn.(^;o-)), (18) 



and AE nsna = - E$\ etc. For the special case when the two initial states ^/ na and 
ty nb are the same and £ = £.', gn a ;n a na reduces to the absolute value of the 2 £ -pole oscillator 
strength 

fL a = {2i T) ^ {Ll] °)\\T,Qi°iYi(*i)\\xn.(L.-, <?))? • (i9) 



Similarly, we have 

B 2 „ C 2 (L,L',L s ,L t ,M) 
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C 2 (L,L',L s ,L t ,M) = ^-G' 2 (L,L > ,L s ,L t ,M)F^(L,L',L s ,L t ), (21) 



with 



G> 2 (L,L>,L s ,L t ,M) = (—l) L+Ls (L, L') 1 / 2 £ ("l) Ms+Mi ^^ 

V-M M s M t )\-M M t MJ 
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Finally, the second-order energy is 

(2) _ ^2n 

R 2n 



vm = _ E ^±)_ (24) 



n>3 

where C 2n (M, ±) are the dispersion coefficients 
C 2n (M,±) = ]T C 1 (L,L',L s ,L t ,M)± E C 2 (L, L' , L s , L t , M) .(25) 

LL'L s Lt LL'L s Lt 
L+L'+2L s +2=2n L+L' +L a +L t +2=2n 

According to third-order perturbation theory, the third-order energy correction is 



v (3) = E E E E E E 
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where Di and D 2 are 



+ E'EE /D (26) 
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D 2 = -(^ )|V|^( ))|(^( )|\/|x„ s (L s M s;< TK t (L,M,;p))| 2 . (28) 
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For the He(2 3 S')-He(2 3 P) system, following a procedure similar to the second-order pertur- 
bation, the third-order energy correction can be expanded in terms of powers of 1/R: 

(3) _ C 9 (M,±) Cll (M,±) 

In this work, we do not present a complete third order calculation. We only consider the 
leading term C 9 (M, ±) of V^ 3 \ which is comparable to the smallest term of V^ 2 \ 

C 9 (M, ±) = C^(M, ±) + G? 2 (M, ±), (30) 

where 



C(?1(Mi±) = E V' v !WW|) x 

y v 5 y ^ ^ ^ (o) (o) 
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< X ^(^;cr)|| ^Q^^C^Oll^Cl;^))^^!^)!! ^ Q^^iCp,) | (p)> , (31) 
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Z7T r 



F D 2 ^ = 37T y, ' ^ngrJX 0, 0, 1, 1) 9nr,n b n b {L t , 1, 1, 1, 1) ^ 
2 n s n t ( A -^n s n a + AE„ tn J 2 |AE„ s „ a AE„ tn J 

The nonrelativistic Hamiltonian for a helium atom in a laboratory frame is 

H = L V ^ - — - — + grage + qn<le + gl (35) 

lulu ^ n 1ui e ^ 1ui e ^ 2 |Rn — Rl| |Rn — R2I | R-l — R-2I 

where m n and m e are the masses of nucleus and electron respectively, q n and q e are their 
charges, and R n , Ri, and R 2 are the corresponding position vectors relative to the laboratory 

6 



frame. Introduce the following transformation 



ri 
r 2 
X 



Rl — Rn , 

R2 — Rn , 
1 



M 7 



m n R„ + m e Ri + m e R 2 ) 



(36) 
(37) 
(38) 



where Mt = m n + 2m e is the total mass of helium atom and X is the position vector for 
the center of mass of helum. In the center-of-mass frame, the Hamiltonian becomes 



H 



1 -Vi-^-Vi-^-V ri .V r , + Mi + M^ + ^ 



(39) 



2/x e 2// e m„ ri r 2 r i2 

where ri 2 = \r\ — r 2 | and fi e = m e m n /(m e + m n ) is the reduced mass between the electron 
and the nucleus. The eigenvalue spectrum and corresponding eigenfunctions are obtained 
by diagonalizing the Hamiltonian ()39|) in the correlated Hylleraas basis set 

{ri ri r\ 2 e~ ar *-»*y™(t u r 2 )} , (40) 

where (f 1, £2) is the coupled spherical harmonics for the two electrons forming a common 
eigenstate of L 2 and L z . Except for some truncations, all terms are included in the basis 
such that 

i + j + k<Q, (41) 

with Q being an integer. As Q increases, the size of basis set is increased progressively. 
It is necessary to transform the transition operator 

T t = q n R e n Y m (Rn) + g e i?%(Ri) + q e RiY eo (R 2 ) (42) 



into the center-of-mass coordinates 



0. 



The transformed Tp with i = 1, 2, and 3 for the 



dipole, quadrupole, and octupole moments for a neutral helium are 

2 
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„m e \ A n /30 m e ,■■>) 
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(43) 
(44) 
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3 /35 
2 V 27 
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2i 



r 2 r 2 ((ri <g> fi) (2) ® f 2 )^ ; + r 2 ri((r 2 ® r 2 ) 



(3) 



»(2) 



(3) 



(45) 
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It is noted that for the case of infinite nuclear mass, the above operators reduce to 



7T = -E^ (r,). (46) 



For the finite nuclear mass case, however, Tg can not be obtained by a simple mass scaling 
from Tf°, except T\ which does not contain m e /Mr explicitly for a neutral system. 

Table gives the convergence pattern of Cq(M,±) for the case of infinite nuclear mass 
as the sizes of basis sets, including the two initial states and the four intermediate states, 
increase progressively. Table ITT1 lists the contributions to Cq(M,±) from three pairs of 
intermediate symmetries ( 3 P, 3 S), ( 3 P,(pp) 3 P) doubly-excited states, and ( 3 P, 3 D). Our 
final results for C 3 (M, ±), C 6 (M, ±), C 8 (M, ±), C 9 (M, ±), and C 10 (M, ±) are presented in 
Table HTT1 To our knowledge, no definitive calculations have been reported for the dispersion 
coefficients C 6 (M, ±), C 8 (M, ±), C 9 (M, ±), and C 10 (M, ±). Table HV| is a comparison with 
the existing values of Cg(M, ±) and C§(M, ±). Our values for C$(M, ±) differ from those 
given in Refs. Q and 5]- The origins of these discrepancies are uncertain. Drake 
tabulated the oscillator strengths of helium, including 3 He and 4 He. In order to extract the 
square of the transition matrix element connecting 2 3 S and 2 3 P, besides some numerical 
constants, the tabulated value should be multipied by 1 + 2m e /m n and divided by the 
transition energy corrected for the finite nuclear mass. Our results are in perfect agreement 
with Drake's values such obtained. For C 6 (M, ±), our values differ from the values used in 
the work of Venturi et al. [5] at the level of 0.8% and 0.9% respectively for Cq(0,±) and 
C R (±1,±). 
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TABLE I: Convergence characteristics of Cq(M, ±), in atomic 
units, for the °°He(2 3 S)-°°He(2 3 P) system. Ns s , N 3p , N* 3s , 
Ns p , N( pp y3 P , and N:3d denote respectively the sizes of bases 
for the two initial states and the four intermediate states of 
symmetries 3 S, 3 P, (pp) 3 P, and 3 D. 

Na s N Sp Nt s m ip N (pp)Sp N, D C 6 (0,±) C 6 (±l,±) 
1330 1360 560 1360 1230 853 2640.233 681 1862.572 368 
1540 1632 680 1632 1430 1071 2640.233 694 1862.572 376 
1771 1938 816 1938 1650 1323 2640.233 700 1862.572 380 



TABLE II: Contributions to Cq(M,±), in atomic units, for 
the °°He(2 3 5)-°°He(2 3 P) system from ( 3 P, 3 S), ( 3 P, {pp) 3 P), 
and ( 3 P, 3 D) symmetries. 



Symmetries C 6 (0,±) C 6 (±l,±) 
( 3 P, 3 S) 684.091969(2) 171.0229919(1) 
( 3 P, (pp) 3 P) 1.31648932(2) 3.29122329(3) 
( 3 P, 3 D) 1954.825248(4) 1688.258168(3) 
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TABLE III: C 3 (M,±), C 6 (M,±), C 8 (M,±), C 9 (Af,±), and 
Ci (M, ±), in atomic units, for the He(2 3 5)-He(2 3 P) system. 



Mass 


^He-^He 


4 He- 4 He 


3 He- 3 He 


(life) +) 


±12.8154931075(4) ±12.8181751205(4) ±12.8190526019(4) 




^6.4077465536(2) 


^6.4090875603(2) 


=F6.4095263011(2) 




2640.2338(1) 


2641.5083(2) 


2641.9255(3) 




1862.5724(1) 


1863.4726(2) 


1863.7674(4) 




311901.2(4) 


311955.4(5) 


311972.8(1) 




lo4lyyo^z ) 


1542352(1) 


1542470(1) 




168906.5(4) 


168921.6(3) 


168926.5(2) 


c 8 (±i,-) 


103017.3(3) 


103039.5(4) 


103046.7(3) 


C 9 (0, ±) 


±512059.227(6) 


±512572.343(6) 


±512740.318(6) 


c 9 (±i,±) 


=Fll7073.536(2) 


^=117199.211(2) 


^117240.354(2) 


Cio(0,+) 


2.922482(3) x 10 7 


2.922304(5) x 10 7 


2.922244(3) x 10 7 


Cio(0,-) 


1.857456(3) x 10 s 


1.8574503(3) x 10 8 1.8574492(3) x 10 8 


Ci (±l,+) 


1.611301(3) x 10 7 


1.611325(3) x 10 7 


1.611334(5) x 10 7 


Ci (±l,-) 


2.40597(3) x 10 6 


2.40608(1) x 10 6 


2.40613(2) x 10 6 



TABLE IV: Comparison of C 3 (M,±) and C 6 (M,±) for the 
4 He(2 3 5)- 4 He(2 3 P) system. 



Author C 3 (0,±) C 3 (±L±) C 6 (0,±) C 6 (±l,±) 

Venturi et at. [5j ±12.82044 =p6. 41022 2620.76 1846.60 

Leonard et al. [2] ±12.810(6) =f6.405(3) 

This work ±12.8181751205(4) =p6. 4090875603(2) 2641.5083(2) 1863.4726(2) 
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